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INTRODUCTION 


In, many wind tunnels, the model support sting drive can both 
pitch and roll the sting and the model. This allows tests to be 
made in both angle of attack, a, and angle of sideslip, p, in the 
wind tunnel. Determination of the sting pitch and roll angles 
necessary to obtain desired a and p angles is fairly simple when 
the sting is straight. If the sting is offset or sting bending 
occurs, however, it can be very difficult to determine the sting 
pitch and roll angles to give the desired o and p on the model. 

In order to solve this problem, a computer program has been 
developed to compute the pitch and roll position of a wind tunnel 
sting to position the model at the desired a and p with respect 
to the tunnel air stream for stings with offset angles in yaw, 
pitch and roll; and sting bending in yaw, pitch and roll. 
Computer programs that have been developed previously have been 
able to calculate the sting position for stings with one or two 
offset angles but generally they have not been able to accept 
sting angles in all three directions and/or stream flow angles. 
Also, previous programs usually generate the cosine of the sting 
pitch or roll angles and very elaborate methods are necessary to 
decide whether the pitch and roll are positive or negative angles 
in order to orientate the model correctly. The computer programs 
that are described in this report generate the sine or tangent of 



angles to reduce the ambiguity of whether the angles are positive 
or negative. 

Two computer programs have been developed and are described 
in this report. These programs cover the case of no 
accelerometers on board the model and the case of accelerometers 
on board the model to measure model pitch and/or model roll with 
respect to gravity. 

These programs are iterative programs in that they calculate 
the a and p of the model based on assumed pitch and roll angles 
of the sting and then calculate new positions for the sting pitch 
and roll to move the model closer to the desired a and p. This 
process is continued until the calculated a and p of the model 
match the desired a and p within a small increment. This 
procedure insures that the a and p calculated by the data 
reduction program based on the available information such as 
stream flow angles, sting bending due to measured model forces, 
and on board measurements of model pitch and roll, if available, 
will match the desired a and p. 

Both of these programs accept three sting offset angles, 
three sting bending angles, and two tunnel flow angles. In 
addition, the second program accepts on board measured pitch 
angle and on board measured roll angle, if available. 
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SYMBOLS 


u Velocity along X axis of model 

v Velocity along Y axis of model 

w Velocity along Z axis of model 

X Longitudinal body axis, positive aft, see figure 1. 

Y Lateral body axis, positive to right, see figure 1. 

Z Vertical body axis, positive upward, see figure 1. 

a Model angle of attack, a = arctan (w/u) , see figure 1. 

P Model angle of slideslip, p = arcsin (-v/Voo) , see 

f igure 1 . 

9 Model pitch angle, positive direction is nose up, see 

figure 1. 

9 of f Sting offset angle in pitch 

9 s k Sting bending in pitch 

9 Sting pitch command 

<t> Model roll angle, positive direction is right wing down, 

see figure 1. 

^off Sting offset angle in roll 

0 sb Sting bending in roll 

<t> Sting roll command 

Model yaw angle, positive direction is nose right, see 
figure 1. 

V» 0 ff Sting offset angle in yaw 

^ sb Sting bending in yaw 
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DETERMINATION OF ANGLE OF ATTACK AND ANGLE OF SIDESLIP 


FOR A WIND TUNNEL MODEL 

The definition of angle of attach, a, on a wind tunnel model 
is given in reference 1 as the arctan (w/u) where w is the 
component of the free stream velocity along the Z axis of the 
model (vertical axis with the positive direction upward) and u is 
the component of the free stream velocity along the X axis of the 
model (longitudinal axis with the positive direction aft) . This 
definition applies no matter what the orientation of the model 
is. By examining the signs of u and w, the correct quadrant for 
a can be determined between -180° and +180°. If both u and w are 
zero, then the angle of attack is indeterminate, but in these 
programs, when both u and w are equal to zero, a is defined to be 
equal to zero. 

The angle of sideslip, p , is defined in reference 1 as the 
arcsin (-v/Voo) where v is the component of the free stream 
velocity along the Y axis of the model (the lateral axis with the 
positive direction out the right or starboard wing) and Voo is the 
total free stream velocity. This means that p is positive when 
the flow is from the right. 

In these computer programs, the free stream velocity, Voo, is 
set equal to one, and the components of the free stream velocity 
along the wind tunnel axis system (u, v and w) are calculated 
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from the angles of upwash and sidewash in the wind tunnel. The 
upwash angle (UWA) is defined to be positive when the flow is 
upward in the tunnel (i.e., the w component is positive) and the 
sidewash to be positive when the flow is from the right to the 
left (from the starboard to the port). For positive wind tunnel 
sidewash angle and the model pitch, roll and yaw angles equal to 
zero, the v component of the model velocity is negative and the 
model sideslip angle, p, is positive. 

The three velocities (u, v and w) along the three axes of 
the model are recalculated for each rotation of the model. By 
convention, the rotations are taken in the order of yaw, pitch 
and then roll (i.e., rotation about the Z axis, the Y axis, and 
then the X axis of the model) . The equations for the u, v and w 
velocities after rotation through each of the angles are given 
below. 

Yaw (rotation about Z axis) , 

UA - U * cos(^) - V * sin(tf) 

VA = V * cos (ip) + U * sin (i>) 

WA - W 

Pitch (rotation about Y axis) , 6 
UB = UA * cos(0) - WA * sin(0) 

VB = VA 

WB - WA * cos ( 0 ) + UA * sin(0) 

Roll (rotation about X axis) , 4 
UC = UB 

VC - VB * cos(*) - WB * sin(^) 

WC - WB * cos(*) + VB * sin(^) 
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Where ip is the angle of yaw (positive for nose right) , 6 is 
the angle of pitch (positive for nose up) , and <f> is the angle of 
roll (positive for right wing down). These angles are shown in 
figure 1 from reference 2. 

The equations given above are used in subroutine ALPBET of 
program STNGOPR to calculate the velocities u, v and w after each 
rotation of the model. After the three velocities are 
determined, the angle of attack, a, and angle of sideslip, p, are 
calculated using the formulas given above. 


DESCRIPTION OF THE COMPUTER PROGRAMS 

Two computer programs were developed to calculate the sting 
pitch and roll angles necessary to position the model to obtain 
the desired a an p in the wind tunnel. 

The first program, STNG, calculates the sting pitch and roll 
based on the sting geometry, such as sting offsets and sting 
bending. This program consists of the main program and two 
subroutines. A listing of the program is given in Appendix A. 

The second program, STNGOPR, calculates the sting pitch and 
roll from the sting geometry but also uses inputs from 
accelerometers on board the model that measure the model pitch 
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and roll angles. This program consists of a main program and 
four subroutines and a listing is given in Appendix B. 

A list and description of the variables used in these 
programs is given after the description of the programs. 


PROGRAM STNG 

Program STNG calculates the sting pitch, 6 sc * and roll, <f> S(i 
necessary to obtain a desired model angle of attack, a, and angle 
of sideslip, p. The program actually first calculates the a and 
p for an assumed sting pitch and roll. Then, by iteration, the 
assumed pitch and roll angles are changed until the calculated a 
and p agree with the desired or command a and p (ALPC and BETC) . 

Initial attempts to simple algorithms to determine which 
direction to move the sting to bring the model to the desired a 
and p, such as changing the sting pitch to change a and changing 
the sting roll to change p did not always converge to the desired 
a and p. Even attempts to control a by the sting movement (pitch 
or roll) that had the most influence on a and then controlling p 
by the remaining sting movement (roll or pitch) were not always 
successful, depending on sting geometry and the relation of the 
starting position of the sting to the final correct position of 
the sting. 
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To overcome these difficulties, two factors which are a 
combination of the errors in a ad ^ were developed for the 
program. These factors, called FT and FP are a linear 
combination of the errors in a and p (DALP and DBET) multiplied 
by a weighting factor which is proportional to the change in a 
and p due to a unit change in the sting pitch command (THESC) and 
the sting roll command (PHISC) respectively. The weighting 
factors for the function FT are FTA and FTB. FTA and FTB are 
proportional to the change in o and the change in p with a unit 
change in sting pitch. The function FT, then, is a linear 
combination of the weighting factor FTA times the error in a 
added to the weighting factor FTB times the error in p. Thus: 

FT = FTA* DALP+FTB* DBET 

After the function FT is found, the program calls subroutine 
CONV which uses two values of sting pitch position (THESC and 
THESCSV) and the two corresponding values of the function FT (FT 
and FTSV) to compute a new value for the sting position which 
should reduce the absolute value of FT. This subroutine uses a 
procedure, which is equivalent to Newton's method of finding the 
roots of an equation, to find the sting position that should 
reduce FT. 

Each time subroutine CONV is called to find the sting 
position to minimize FT, subroutine VEL is called to determine 
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a 




the a and p at the new sting pitch angle. These new values for a 
and p are then used to calculate a new FT. This process is 
repeated until the absolute value of FT is less than .00001 
(TOLF) or for a maximum of three iterations. 

After FT is minimized, the weighting factor for the sting 
roll function (FP) are found. These weighting factors (FPA and 
FPB) are proportional to the change in a and the change in p, 
respectively, with a change in sting roll command (PHISC) . The 
roll function, FP, then is: 

FP = FPA* DALP+ FPB* DBET 

In the same manner as with FT, subroutine CONV is called to 
find a sting roll command (PHISC) that reduces the value of FP. 

After FT and FP have been successively reduced, the entire 
procedure is repeated until the sum of the errors in a 
(ABS(DALP)) and p (ABS(DBET)) is less than the tolerance (TOLDAB) 
or for a maximum of six iterations. 

In general, there are two positions of the sting drive that 
satisfy the command a and p. There are certain special cases 
where there are infinitely many positions and there are certain 
other cases where the sting cannot be positioned by pitch and 
roll drives to satisfy the commanded a and p. In those cases 
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where there are two positions, they are nearly at a 180° rotation 
around the X axis of the wind tunnel to each other such that in 
one position the model is upright and in the other the model is 
inverted. The solution developed by the computer program can be 
either the upright or the inverted solutions depending upon the 
initial position of the sting and the a and p desired. The 
program checks the sting roll position of the solution to 
determine if it is outside the limits of -85° to 100°, and if so 
the program is rerun with the initial conditions of pitch equal 
to the negative of first solution pitch, and roll equal to the 
first solution roll +180°. The solution with the new initial 
values for sting pitch and sting roll is very quick and the sting 
roll is within the limits of -85° to +100°. If an inverted 
solution is desired, the limits on the sting roll position can be 
changed to +80° for the lower limit and +280° for the upper 
limit. 

The following is a list and description of the variables 
used in program STNG (all angles are degrees) : 

ALP Angle of attack of model, a. 

ALPC Command or desired angle of attack. 

BET Angle of sideslip model, p. 
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BETC 


Command or desired angle of sideslip. 


DALP 

DBET 

FP 

FPA 

FPB 

FT 

FTA 

FTB 

PHIOFF 

PHISB 


Difference between angle of attack of the model and 
desired angle of attack. 

Difference between angle of sideslip of the model and 
desired angle of sideslip. 

Function of sting roll ( 4 SC ) which is to be minimized by 
rolling the sting. 

Weighting factor for DALP in the function FP. 

Weighting factor for DBET in the function FP. 

Function of sting pitch (0 which is to be minimized 
by pitching the sting. 

Weighting factor for DALP in the function FT. 

Weighting factor for DBET in the function FT. 

Sting offset angle in roll 
Sting bending in roll (tf sb ) 
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PHISC 


PHISLL 


PHISLU 


PSIOFF 


PSISB 


SWA 


THEOFF 


THESB 


THESC 


TOLDAB 


TOLF 


Sting roll command (0 ) required to position the model 

SO 

at the desired a and p. 

Lower limit on sting roll. 

Upper limit on sting roll. 

Sting offset angle in yaw • 

Sting bending in yaw (^g^) • 

Wind tunnel free stream sidewash angle, positive for 
flow from right. 

Sting offset angle in pith (# Q f f ) • 

Sting bending angle in pitch ( 6 . 

Sting pitch command ( 6 ^ required to position the model 
at the desired a and p. 

Convergence tolerance for the sum of the absolute values 
of DALP and DBET. 

Convergence tolerance for the value of FT and FP. 
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UWA Wind tunnel free stream upwash angle, positive for flow 

from below. 

u Free stream velocity component in the longitudinal 

direction in the wind tunnel (the total free stream 
velocity is assumed to be 1.0). 

v Free stream velocity component in the lateral direction 

(flow from the right when looking forward is positive). 

w Free stream velocity component in the vertical direction 

(upward flow is positive) . 


SUBROUTINE VEL 

The purpose of subroutine VEL is to calculate the angle of 
attack, a, and angle of sideslip, p, of the wind tunnel model. 
The subroutine requires inputs of the velocity components in the 
wind tunnel u, v and w; the sting drive angles; the sting offset 
angles; and the sting bending angles. The subroutine calculates 
the components of the free stream velocity along the three axes 
of the model after each rotation angle using the formulas given 
in the section "Determination of Angle of Attack and Angle of 
sideslip for a Wind Tunnel Model." After the last rotation, 
these velocities are used to calculate the angle of attack, a, 
and angle of sideslip, p, using the following formulas: 
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a = arctan (w/u) 


p = arcs in (-v/Voo) 

Voo is set equal to one in the main program (STNG) and therefore: 

p = arcsin (-v) 

The following is a list and description of the additional 
variables used in subroutine VEL: 


UB, . . . ,UI 


Longitudinal velocity in the model axis 
system after each rotation. 


VB, . . . ,VI 


Lateral velocity in the model axis system 
after each rotation. 


m, ... ,wi 


Vertical velocity in the model axis system 
after each rotation. 
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SUBROUTINE CONY 


The purpose of subroutine CONV is to minimize a function Y = 
F(X). The function is calculated in the calling program and the 
subroutine calculates a new value of the independent variable X 
that will make the value of the dependent variable Y, nearer to 
zero. The new value of X is determined by calculating where a 
straight line through the last two previous pairs of points; X 
and Y, and XSAVE and YSAVE ; intersects the X axis. 

This new value is returned to the calling program in the X 
parameter and the old X and Y are placed in the XSAVE and YSAVE 
parameters at the end of the subroutine. In order to improve the 
stability of this procedure, certain limits are placed on the 
distance the new X value can be from the old X and XSAVE values. 
Line 14 of the subroutine CONV limits the new X value to be no 
more than 3.5 times the distance between the old X and XSAVE 
values away from the average of these values. Also if the two 
previous values for X are the same, the new X is calculated in a 
special way in statement 3. If the two previous Y values are the 
same, the new X is calculated to be the average of the two 
previous X's in statement 4. 
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The following is a list and description of the variables used 

in subroutine CONV: 

DX Absolute distance between XI and X2. 

X Latest value of the independent variable. Also returned 
as the next value to be tried for the independent 
variable. 

XA Average of XI and X2 . 

XK Slope of line between XI, Y1 and X2, Y2. 

XSAVE Previous value of independent variable. Also returned 
as previous value of X. 

XI Previous value of the independent variable (same as 
XSAVE) . 

X2 Latest value of the independent variable. 

Y Latest value of the dependent variable. 

YSAVE Previous value of the dependent variable. Also returned 
as the previous value of Y. 
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Y1 Previous value of the dependent variable (same of 

YSAVE) . 

Y2 Latest value of the dependent variable. 

PROGRAM STMGOPR 

Program STNGOPR calculates the sting pitch and roll angles 
reguired to obtain the desired angle of attack and angle of 
sideslip in the wind tunnel when the model has on board 
measurements of the model pitch and roll (or pitch only if the 
roll is not available) . The model pitch and roll is assumed to 
be measured relative to gravity. 

At first it would appear that these cases would simplify the 
calculation of the alpha and beta of the model in the wind tunnel 
since the effects of sting offsets and sting bending are already 
included in the on board measurements of model pitch and roll. 
This is true in the case where a straight sting is used and the 
model yaw is assumed to equal zero. Since most wind tunnel tests 
use straight stings, on board measurements of model pitch and 
roll is a very valuable method of determining model alpha and 
beta in the wind tunnel. For bent stings, however, model yaw 
cannot be assumed to equal zero and model pitch and roll alone is 
not sufficient to determine model angle of attack and angle of 
sideslip. 
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In order to make program STNGOPR more generally applicable, 
no assumptions were made about the sting offset angles or the 
model yaw angle in the wind tunnel. Therefore, the program is 
applicable not only to the case of straight sting, but also to 
the case where the sting is offset and/or where sting bending 
occurs . 

Program STNGOPR is similar to program STNG in that they both 
calculate the model a and 0 at a given sting position and then, 
by the use of a factor, try to reduce the difference between the 
calculated a and 0. 

Program STNGOPR uses the same functions (FT and FP) as 
program STNG uses to reduce the error between the calculated a 
and 0 and the desired a and 0. The major difference between the 
two programs is that STNGOPR uses a three-step process to 
calculate a and 0 instead of the one step that STNG uses 
(subroutine VEL) . The three steps are: first, calculate the 

yaw, pitch and roll angles of the model from the sting support 
system angles and the sting geometry (subroutine SIMUST) , second, 
correct the pitch and roll of the model by the difference between 
the measured model pitch and roll, and the calculated model pitch 
and roll ( DTHEMOB and DPHIMOB) determined at the beginning of the 
program, and third, calculate the a and 0 of the model from the 
model yaw, pitch and roll angles (subroutine ALPBET) . 
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When programs STNG and STNGOPR are used in actual wind tunnel 
situations, the sting drive angles calculated by these programs 
will change somewhat as the sting and model move to the commanded 
positions because the sting bending angles will change as the 
angles of attack and sideslip of the model change and because the 
differences between the calculated model pitch and roll and the 
on board measured model pitch and roll change as the model 
attitude changes. The final position, however, will be the 
correct position to obtain the desired a and p, since the final 
values for the sting bending and measured model pitch and roll 
will be the same as those used by the wind tunnel data reduction 
program . 

The following is a list and description of the additional 
variables used in program STNGOPR (all angles are in degrees) : 

DPHIMOB Difference between PHIMOB and PHIMT 

determined at the beginning of the 
program. 

DTHEMOB Difference between THEMOB and THEMT 

determined at the beginning of the 
program . 

PHIMOB Roll of model as measured by on board 

accelerometers . 
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PHIMOBT 


Theoretical on board roll of model. 


PHIMOBT = PHIMT + DPHIMOB 


PHIMT 


PHIS 


PSIMT 


THEMOB 


Theoretical roll of model as determined by 
subroutine SIMUST from the sting drive 
angles, sting offset angles, and sting 
bending angles. 

Actual sting drive roll angle. 

Theoretical yaw angle of the model as 
determined by subroutine SIMUST from the 
sting drive angles, sting offset angles, 
and sting bending angles. 

Pitch of model as measured by on board 
accelerometers . 


THEMOBT 


Theoretical on board pitch of model. 
THEMOBT = THEMT + DTHEMOB 
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THEMT 


Theoretical pitch of model as determined 
by subroutine SIMUST from the sting drive 
angles, sting offset angles and sting 
bending angles. 

THES Actual sting drive pitch angle. 

SUBROUTINE SIMUST 

Subroutine SIMUST simulates the sting-support-sting-model 
system mathematically to calculate the model yaw, pitch and roll 
angles from the sting drive, sting offset and sting bending 
angles. The method of calculating the model yaw, pitch and roll 
is discussed in references 2 and 3, and is explained below. 

The pitch angle of the model is determined by calculating the 
X component, in the model axis system, of a unit vector in the Z 
direction of the tunnel axis system (XZ) . The pitch angle is 
then the arcsin (-XZ) . The pitch angle can range from -90° to 
+90°. 

The roll of the model is determined by calculating the Y and 
Z components, in the model axis system, of a unit vector in the Z 
direction of the tunnel axis system (YZ and ZZ) . The roll of the 
model is then the arctan (-YZ/ZZ) where the guadrant of the roll 
angle is determined by the signs of YZ and ZZ individually. The 
roll angle can range from -180° to 180°. If both YZ and ZZ are 
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zero (i.e., the pitch angle is ±90°) then the roll of the model 
is determined by the arctan (-YX, ZX) and the yaw of the model is 
defined to be equal to zero (where YX ad ZX are the Y component 
and the Z component respectively of a unit vector in the X 
direction of the tunnel axis system) . 

The yaw of the model is determined by calculating the X 
component in the model axis system of a unit Y vector in the wind 
tunnel axis system, XY, and the X component in the model axis 
system of a unit X vector in the tunnel axis system, XX. The yaw 
of the model is then the arctan (-XY/XX) and can range from -180° 
to 180°. 

The following is a list and description of the additional 
variables used in subroutine SIMUST: 

PHIS Sting drive roll angle 

THES Sting drive pitch angle 

XX, YX, ZX The X, Y and Z components in the model 

axis system of a unit vector in the X 
direction in the tunnel axis system. 
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XY, YY, ZY 


The X , Y and Z components in the model 
axis system of a unit vector in the Y 
direction in the tunnel axis system. 

XZ , YZ , ZZ The X, Y and Z components in the model 

axis system of a unit vector in the Z 
direction in the tunnel axis system. 


SUBROUTINE COMP 


Subroutine COMP is used to calculate the components, in the 
model axis system, of a unit vector in the wind tunnel axis 
system so that the model yaw, pitch and roll angles can be 
determined. This subroutine is very similar to subroutine VEL in 
program STNG which was used to calculate the velocity components 
in the model body axis system. Although this subroutine can 
calculate the X, Y and Z components, in the model axis system, of 
an arbitrary vector in the wind tunnel axis system, it is only 
used in this program to calculate the components of a unit vector 
in the X, Y or Z direction in the wind tunnel axis system. This 
means that one of the components, X, Y or Z, is set equal to one 
and the other components are set equal to zero in the calling 
program argument list. The components in the model body axis 
system of the unit vectors are then used in subroutine SIMUST to 
calculate the model yaw, pitch and roll. 
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The following is a list and description of the additional 


variables used in subroutine COMP: 


X, Y, Z 


XB, . . . ,XI 


YB, . . . , YI 


ZB, • • • , ZI 


X, Y and Z components of a vector in the 
tunnel axis system. 

X component in the model axis system of a 
vector in the tunnel axis system after 
each rotation. 

Y component in the model axis system of a 
vector in the tunnel axis system after 
each rotation. 

Z component in the model axis system of a 
vector in tunnel axis system after each 
rotation. 
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SUBROUTINE ALPBET 


The purpose of subroutine ALPBET is to calculate the angle of 
attack, a, and angle of sideslip, 0, of the wind tunnel model. 

The subroutine requires inputs of the velocity components in the 
wind tunnel, U, V and W, and the model Euler angles, yaw (^») , 
pitch (0) , and roll (<f >) . The subroutine calculates the 
components of the free stream velocity along the three axes of 
the model after each rotation angle using the formulas given in 
the section "Determination of Angle of Attack and Angle of 
Sideslip for a wind Tunnel Model." After rotation through the 
three Euler angles, the velocities are used to calculate the 
angle of attack, a, and angle of sideslip, 0, using the following 
formulas: 


a = arctan (w/u) 

0 = arcs in (-v/Voo) 

Voo is set to one in the main program (STNGOPR) and therefore: 

0 - arcsin (-v) 

The following is a list and description of the additional 
variables used in subroutine ALPBET (all angles are in degrees) : 
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ALP 

BET 

PHIM 

PSIM 

THEM 

UA, UB, UC 
VA, VB, VC 
WA, WB, WC 


Angle of attack of model, a. 

Angle of sideslip of model, p. 

Angle of roll of model, <l>. 

Angle of yaw of model, 0. 

Angle of pitch of model, 9. 

Longitudinal velocity component in the 
model axis system after each Euler angle 
rotation. 

Lateral velocity component in the model 
axis system after each Euler angle 
rotaton. 

Vertical velocity component in the model 
axis system after each Euler angle 
rotation. 
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CONCLUDING REMARKS 


Two programs have been developed to calculate the pitch and 
roll angles of a wind tunnel sting drive system that will 
position a model at the desired angle of attack and angle of 
sideslip in the wind tunnel. These programs account for the 
effects of sting offset angles, sting bending angles and wind 
tunnel stream flow angles. In addition, the second program 
incorporates inputs from on board accelerometers that measure 
model pitch and roll with respect to gravity. 

These program solve for the desired sting pitch and roll with 
an iterative procedure using the forward equations that calculate 
the model a and p from the sting geometry and the sting pitch and 
roll. This procedure avoids the ambiguity that is found in many 
inverse solutions that solve for the sting pitch and roll from 
the model a and p. Also, more sting offset angles, sting bending 
angles, and stream flow angles can be taken into account by using 
the forward equations. 

A copy of the source code of these two programs can be 
obtained from the Langley Computer Center with the following 
statements : 


GET, STNG/UN - 69 02 5 ON 
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or GET, STNGOPR/UN = 690250N 


The run times for these programs vary depending upon the 
number of iterations reguired to converge to a solution. When 
compiled and run under Fortran 5 (Fortran 77) on the Control Data 
Corporation Cyber CY180-860 computer at Langley Research Center, 
the run times for STNG is from 0.005 seconds for one iteration to 
0.100 seconds for 36 iterations (the maximum allowed in these 
programs). The run times for STNGOPR range from 0.027 seconds to 
0.352 seconds for one to 36 iterations respectively. 
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APPENDIX A 


COMPUTER LISTING OF PROGRAM STNG 


This appendix contains a computer listing of the program STNG 
which calculates the wind tunnel sting pitch and roll angles 
required to obtain the angle of attack, a, and angle of sideslip, 
p, on a wind tunnel model. The program accepts stream flow 
angles in two directions and sting offsets and sting bending in 
three directions. 
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PROGRAM STNG (INPUT, OUTPUT) 

C THIS PROGRAM CALCULATES THE STING DRIVE PITCH AND ROLL ANGLES TO 

C POSITION A WIND TUNNEL MODEL AT THE COMMANDED ANGLE OF ATTACK (ALPC) 

C AND ANGLE OF SIDESLIP (BETC). THE PROGRAM ACCOUNTS FOR STING OFFSET 

C ANGLES, STING BENDING ANGLES AND STREAM FLOW ANGLES. 

C 

C CODED BY — JOHN B. PETERSON, JR. NASA/LARC/TAD/NTFOB 1988 

C STING OFFSETS 

PSIOFF=0. 

THEOFF-45. 

PHIOFF-20. 

C STREAM FLOW ANGLES 

C POSITIVE FOR FLOW FROM BELOW AND FROM RIGHT 

SWA*. 0 
UWA* . 1 

C FREE STREAM VELOCITIES (TOTAL VEL. - 1.0) 

U=SQRT( 1./ (l.+(TAND(SWA))**2+(TAND(UWA))**2 ) ) 

V=-U*TAND(SWA) 

W= U*TAND(UWA) 

C INITIAL VALUES 

THESC=0.0 
PHISC*0.0 

C LIMITS 

PHISLL=-85 . 

PHISLU=100. 

FOR INVERTED RUNS USE FOLLOWING LIMITS 
PHISLL-80. 

PHISLU*265 . 

TOLF=. 00001 
TO LDAB=. 00001 

C INPUTS TO CONTROL PROGRAM 

C COMMAND ANGLES 

ALPC* 5. 

BETC- 10. 

C STING DEFLECTIONS 

PSISB-0.4 
THESB-0.4 
PHISB-0.2 


C 


END OF INPUTS TO CONTROL PROGRAM 


50 CONTINUE 


PRINT 99 

PRINT 98 , THESC , PHI SC , ALPC , BETC , PS IOFF , THEOFF , PHIOF F 

C CONVERGE ON ALPC AND BETC 

DO 100 ICONV=l ,6 
PRINT 96 
PRINT 97 

THESCSV=THESC-1. 

CALL VEL (U,V,W, PS IOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESCSV ,PHISC ,ALPSV,BETSV) 

DALPSV=ALPSV-ALPC 

DBETSV=BETSV-BETC 

CALL VEL ( U , V , W , PSIOFF , THEOFF , PHIOFF , PSISB ,THESB , PHISB , 

* THESC, PHISC, ALP, BET) 

DALP=ALP-ALPC 

DBET=BET-BETC 

FTA=ALP -ALPSV 
FTB=BET -BETSV 
X=1 . / (ABS(FTA)+ABS(FTB) ) 

FTA°FTA*X 

FTB=FTB*X 

h'TSV = FTA*DALPSV+FTB*DBETSV 
FT =FTA*DALP +FTB*DBET 
DO 200 ICTHE=1 , 3 

CALL C0NV( THESC, FT, THESCSV, FTSV) 

CALL VEL (U,V,W, PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESC, PHISC, ALP, BET) 

DALP=ALP-ALPC 

DBET=BET-BETC 

FT=FTA*DALP+FTB*DBET 

PRINT 98 , THESC , PHISC , DALP , DBET , FT 
IF(ABS(FT) .LT.TOLF) GO TO 210 
200 CONTINUE 
210 CONTINUE 

I F ( A BS ( DALP )+ABS ( DBET ) . LT . TOLDAB ) GO TO 1000 
IF (THESC. EQ.O. ) THESC=. 000001 
PRINT 96 
PRINT 97 
PHISCSV=PHISC-1. 

CALL VEL (U,V,W, PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESC, PHISCSV, ALPSV, BETSV) 

DALPSV=ALPSV-ALPC 

DBETSV=BETSV-BETC 

CALL VEL (U,V,W, PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESC, PHISC, ALP, BET) 

DAL P=AL P-ALPC 
DBET=BET-BETC 
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FPA=ALP-ALPSV 
FPB^BET-BETSV 
DENOM-ABS (FPA )+ABS ( F PB ) 

IF(DENOM.EQ.O. ) GO TO 310 

X=l . /DENOM 

FPA“FPA*X 

FPB-FPB*X 

FPSV*FPA*DALPSV+FPB*DBETSV 

FP -FPA*DALP +FPB*DBET 

DO 300 ICPHI*1 , 3 

CALL CONV ( PHI SC , FP , PHISCSV , FPSV ) 

CALL VEL (U, V,W ,PSIOFF ,THEOFF,PHIOFF ,PSISB ,THESB , PHISB , 
* THESC , PHISC , ALP , BET ) 

DALP=ALP-ALPC 

DBET-BET-BETC 

FP=FPA*DALP+FPB*DBET 

PRINT 98, THESC, PHISC, DALP,DBET,FP 
IF(ABS(FP) . LT.TOLF) GO TO 310 
300 CONTINUE 
310 CONTINUE 

IF(ABS(DALP)+ABS(DBET) .LT.TOLDAB) GO TO 1000 
100 CONTINUE 
1000 CONTINUE 

C CHECK TO SEE IF CONVERGED OUTSIDE PHIS LIMIT 
IF( PHISC. GT. PHI SLU) THEN 
PHISC=PHISC-180. 

THESC— THESC 
GO TO 50 
END IF 

I F( PHISC.LT. PHIS LL) THEN 
PHISC=PHISC+180. 

THESC— THESC 
GO TO 50 
END IF 

PRINT 96 
PRINT 95 

PRINT 98, THESC, PHISC, ALP, BET 
STOP 


99 

FORMAT (" 

THESC 

PHISC 

ALPC 

BETC 

PSIOFF 


* 

THEOFF 

PHIOFF") 




98 

FORMAT ( 

8F10.5) 





97 

FORMAT (" 

THESC 

PHISC 

DALP 

DBET 

F") 

96 

FORMAT ( ) 






95 

FORMAT ( " 

THESC 

PHISC 

ALP 

BET”) 



END 



SUBROUTINE VEL (U, V,W,PSI0FF,THE0FF,PHI0FF,PSIS8,THESB,PHISB, 
* THESC ,PHISC ,ALP,BET) 

DOR-57.2957795 

C STING PITCH (Y) 

UB-U*COSD(THESC)-W*SIND(THESC) 

VB-V 

WB=W*COSD( THESC )+U*SIND( THESC) 

C STING ROLL (X) 

UC-UB 

VC-VB*COSD(PHISC)-WB*SIND(PHISC) 

WC-WB*COSD(PHISC)+VB*SIND(PHISC) 

C OFFSET YAW (Z) 

UD-UC*COSD(PSIOFF)-VC*SIND(PSIOFF) 
VD=VC*COSD(PSIOFF)+UC*SIND(PSIOFF ) 

WD-WC 

C OFFSET PITCH (Y) 

UE=UD*COSD(THEOFF)-WD*SIND(THEOFF) 

VE-VD 

WE»WD*COSD(THEOFF )+UD*SIND(THEOFF) 

C OFFSET ROLL (X) 

UF-UE 

VF-VE*COSD(PHIOFF)-WE*SIND(PHIOFF) 

WF-WE*COSD(PHIOFF)+VE*SIND(PHIOFF) 

C STING BENDING IN YAW (Z) 

UG-UF*COSD(PSISB)-VF*SIND(PSISB) 
VG»VF*COSD(PSISB)+UF*SIND(PSISB) 

WG-WF 

C STING BENDING IN PITCH (Y) 

UH-UG*COSD(THESB)-WG*SIND(THESB) 

VH-VG 

WH=WG*COSD(THESB)+UG*SIND(THESB) 

C STING BENDING IN ROLL (X) 

UI-UH 

VI=VH*COSD(PHISB)-WH*SIND(PHISB) 

WI=WH*COSD(PHISB)+VH*SIND(PHISB) 

C ALPHA AND BETA 

IF(WI . EQ. 0. . AND.UI .EQ.O. )UI=. 0000001 
ALP=ATAN2(WI ,UI)*DOR 
IF (VI .LT .-1 . )VI— 1. 

IF(VI.GT.l.) VI-1. 

BET-ASIN(-VI)*DOR 

RETURN 

END 
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SUBROUTINE CONV (X, Y.XSAVE.YSAVE) 

XI=XSAVE 

Y1=YSAVE 

X2=*X 

Y2=Y 

IF(X2 . EQ.X1 ) GO TO 3 
I P(Y2 . EQ. Y1 ) GO TO 4 
XK-(Y2-Y1)/(X2-XL) 

X-X2-Y2/XK 

XA-(Xl+X2)/2. 

DX=ABS(X2-X1 ) 

IF(ABS(X-XA).GT.3.5*DX) X=XA+3 . 5*SIGN( DX, (X-XA) ) 
GO TO IOO 

3 X*X2-Y2 
GO TO 100 

4 X-(X2+Xl)/2. 

100 XSAVE-X2 

YSAVE-Y2 

RETURN 

END 


35 



APPENDIX B 


COMPUTER LISTING OF PROGRAM STNGOPR 


This appendix contains a computer listing of program STNGOPR 
which calculates the wind tunnel sting pitch and roll angles 
required to obtain the angle of attack, a, and angle of sideslip, 
p, on a wind tunnel model with on board accelerometers to measure 
model pitch and/or roll angles. The program accepts the 
accelerometer measurements of model pitch and roll, wind tunnel 
stream flow angles in two directions, and sting offsets and sting 
bending in three directions. 
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PROGRAM STNGOPR(INPUT, OUTPUT) 

C THIS PROGRAM CALCULATES THE STING DRIVE PITCH AND ROLL ANGLES TO 
C POSITION A WIND TUNNEL MODEL AT THE COMMANDED ANGLE OF ATTACK (ALPC) 

C AND ANGLE OF SIDESLIP (BETC). THE PROGRAM ACCEPTS INPUTS FROM ACCELEROMETERS 
C ON BOARD THE MODEL TO MEASURE THE MODEL PITCH AND ROLL RELATIVE TO GRAVITY 
C AND IT ACCOUNTS FOR STING OFFSET ANGLES, STING BENDING AND STREAM FLOW ANGLES. 
C 

C CODED BY ~ JOHN B. PETERSON, JR. NASA/LARC/TAD/NTFOB 1988 

C STING OFFSETS 

PSI0FF=0. 

THEOFF=45. 

PHI0FF=20. 

C STREAM FLOW ANGLES 

C POSITIVE FOR FLOW FROM BELOW AND FROM RIGHT 

SWA-.O 
UWA= . 1 

C FREE STREAM VELOCITIES (TOTAL VEL. =1.0) 

U=SQRT( 1./ ( 1 .+(TAND(SWA) )**2+(TAND(UWA) )**2 ) ) 

V=-U*TAND(SWA) 

W= U*TAND(UWA) 

C INITIAL VALUES 

THESC=0.0 
PHISC=0.0 
THEM0B=0.0 
PHIM0B=0.0 

LIMITS 
PHISLL=-85. 

PHISLU=100. 

INVERTED RUNS USE FOLLOWING LIMITS 
PHISLL=80. 

PHISLU=265. 

TOLF=. 00001 
TOLDAB=. 00001 

50 CONTINUE 

C INPUTS TO CONTROL PROGRAM 

C COMMAND ANGLES 

ALPC=5. 

BETC=10. 

C STING DEFLECTIONS 

PSISB=0.4 


C 

FOR 

C 
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THESB=0.4 

PHISB=0.2 

C STING POSITION 

THES=THESC 

PHIS=PHISC > 

C MODEL PITCH AND ROLL 

CALL S IMUST( PSIOFF .THEOFF .PHIOFF , PSISB , THESB , PHISB , 
* THES , PHIS ,PSIMOB,THEMOB .PHIMOB) 

THEMO B =THEM0B+0 . 

PHIM0B=PHIM0B+0. 

C END OF INPUTS TO CONTROL PROGRAM 


PRINT 99 

PRINT 98, THESC, PHISC, ALPC.BETC .PSIOFF .THEOFF, PHIOFF 

C START OF CONTROL PROGRAM 

C DETERMINE THEORETICAL PITCH AND ROLL OF MODEL 

C AND COMPARE WITH MEASURED PITCH AND ROLL TO GET ERROR 

CALL SIMUST(PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THES , PHIS , PSIMT ,THEMT , PHIMT)- 
DTHEMOB-THEMOB-THEMT 
DPHIMOB-PHIMOB-PHIMT 

************************************************************ 

C IF THEMOB OR PHIMOB IS NOT AVAILABLE, DTHEMOB OR DPHIMOB 
C SHOULD BE SET TO ZERO 
DTHEMOB-O. 

DPHIMOB-O. 

************************************************************ 

C CONVERGE ON ALPC AND BETC 

PRINT 96 
PRINT 94 

PRINT 98 , THEMO B , PHIMOB , DTHEMOB , DPHIMOB 

DO 100 ICONV=l , 6 
PRINT 96 
PRINT 97 

THKSCSV-THESC-1. 

CALL SIMUST( PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESCSV , PHISC , PSIMT ,THEMT ,PHIMT) 

THEMO BT-THEMT+DTHEMOB 
PHIMOBT=PHIMT+DPHIMOB 

CALL ALPBET(U , V , W , PSIMT , THEMOBT , PHIMOBT , ALPSV , BETSV ) 
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DALPSV-ALPSV-ALPC 

DBETSV-BETSV-BETC 

CALL S IMUST( PS IOF F , THEOFF , PHIOFF , PS ISB,THESB, PHISB, 

* THESC , PH I SC, PS IMT, THEMT, PHIMT) 
THEMOBT=THEMT+DTHEMOB 
PHIMOBT-PHIMT+DPHIMOB 

CALL ALPBET ( U , V , W , PS IMT , THEMOBT , PHIMOBT , ALP , BET ) 

DALP=ALP-ALPC 

DBET=BET-BETC 

FTA=ALP -ALPSV 

FTB=BET -BETSV 

X=1./(ABS(FTA)+ABS(FTB)) 

FTA=FTA*X 

FTB=FTB*X 

FTSV=FTA*DALPSV+FTB*DBETSV 

FT =FTA*DALP +FTB*DBET 

DO 200 ICTHE=1 ,3 

CALL CONV ( THESC , FT , THESCS V , FTSV ) 

CALL SIMUST( PSIOFF ,THEOFF ,PHIOFF ,PSISB , THESB ,PHXSB , 

* THESC , PHISC, PS IMT, THEMT, PHIMT) 
THEMOBT=THEMT+DTHEMOB 
PHIMOBT=PHJMT+DPHIMOB 

CALL ALPBET ( U , V , W , PS IMT , THEMOBT , PHIMOBT , ALP , BET ) 

DALP=ALP-ALPC 

DBET=BET-BKTC 

FT=FTA*DALP+FTB*DBET 

PRINT 98 , THESC , PHISC , DALP , DBET , FT 
IF(ABS(FT) .LT.TOLF) GO TO 210 
200 CONTINUE 
2LO CONTINUE 

I F ( ABS ( DALP )+ABS ( DBET ) . LT . TOLDAB ) GO TO 1000 
IF(THESC.EQ.O. ) THESC=. 000001 
PRINT 96 
PRINT 97 

PHISCSV-PHISC-1. 

CALL SIMUST( PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESC , PHISCSV , PSIMT , THEMT , PHIMT ) 
THEMOBT-THEMT+DTHEMOB 
PHIMOBT-PHIMT+DPHIMOB 

CALL ALPBET ( U , V , W , PS IMT , THEMOBT , PHIMOBT , ALPSV , BETSV ) 

DALPSV-ALPSV-ALPC 

DBETSV=BETSV-BETC 

CALL SIMUST( PSIOFF, THEOFF, PHIOFF, PSISB, THESB, PHISB, 

* THESC , PHISC, PSIMT, THEMT, PHIMT) 
rHEMOBT*THEMT+DTHEMOB 
PHIMOBT=PHIMT+DPHIMOB 

CALL ALP BET ( U , V , W , '*S IMT , THEMOBT , PHIMOBT , ALP , BET ) 

DALP-ALP-ALPC 

DBET=BET-BETC 

FPA=ALP-ALPSV 



FPB-BET-BETSV 
DENOM=ABS ( FPA )+ABS (FPB) 

IF(DENOM.EQ.O. ) GO TO 310 

X=1 . /DENOM 

FPA=FPA*X 

FPB=FPB*X 

FPSV=FPA*DALPSV+FPB*DBETSV 
FP =FPA*DALP +FPB*DBET 
DO 300 ICPHI=1 , 3 

CALL CONV ( PHI SC , FP , PHISCSV , FPSV ) 

CALL SIMUST(PSIOFF ,THEOFF ,PHIOFF ,PSISB ,THESB ,PHISB , 
* THESC , PHISC , PSIMT , THEMT , PHIMT ) 

THEMOBT=THEMT+DTHEMOB 
PHIMOBT=PHIMT+DPHIMOB 

CALL ALPBET ( U , V , W , PSIMT , THEMOBT , PHIMOBT , ALP , BET ) 

DALP=ALP-ALPC 

DBET=BET-BETC 

FP=FPA*DALP+FPB*DBET 

PRINT 98 , THESC , PHISC , DALP , DBET , FP 
IF(ABS(FP) .LT.TOLF) GO TO 310 
300 CONTINUE 
310 CONTINUE 

I F ( ABS (DALP )+ABS ( DBET ) . LT . TOLDAB ) GO TO 1000 
100 CONTINUE 
1000 CONTINUE 

C CHECK TO SEE IF CONVERGED OUTSIDE PHISC LIMIT 
IF(PHISC.GT.PHISLU) THEN 
PHISC=PHISC-180. 

THESC— THESC 
GO TO 50 
END IF 

IF(PHISC.LT.PHISLL) THEN 
PHISC=PHISC+180. 

THESC— THESC 
GO TO 50 
END IF 

PRINT 96 
PRINT 95 



PRINT 98 
STOP 

, THESC, 

PHISC, ALP, 

BET 



99 

98 

FORMAT ( " 

* 

FORMAT ( 

THESC 

THEOFF 

8F10.5) 

PHISC 

PHIOFF”) 

ALPC 

BETC 

PSIOFF" , 

97 

96 

FORMAT ( " 
FORMAT ( ) 

THESC 

PHISC 

DALP 

DBET 

F") 

95 

FORMAT (" 

THESC 

PHISC 

ALP 

BET") 


94 

FORMAT ( " 
END 

THEMOB 

PH I MOB 

DTHEMOB 

DPHIMOB" ) 
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SUBROUTINE SIMUST (PSIOFF,THEOFF,PHIOFF,PSISB,THESB, PHISB, 

* THES, PHIS, PS IMT, THEM!, PHIMT) 

C THIS SUBROUTINE SIMULATES THE STING. IT CALCULATES THE THEORETICAL MODEL 
C YAW, PITCH AND ROLL GIVEN INPUTS OF THE 
C STING OFFSETS, STING BENDING, AND STING PITCH AND ROLL. 

DOR-57.2957795 

CALL COMP (I. ,0. ,0. ,PS10FF,THE0FF,PHI0FF,PSISB,THESB, PHISB, 

* THES, PHIS, XX, YX,ZX) 

CALL COMP (0 . , 1 . ,0. , PSIOFF ,THEOFF ,PHIOFF ,PSISB ,THESB , PHISB , 

* THES, PHIS, XY,YY,ZY) 

CALL COMP (0. ,0. ,1. ,PSIOFF,THEOFF,PHIOFF,PSISB,THESB, PHISB, 

* THES, PHIS, XZ,YZ,ZZ) 

C CALCULATE PITCH OF MODEL 

IF(XZ.LT.-1. )XZ=-1. 

IF(XZ.GT. l.)XZ= 1. 

THEMT-AS IN ( -XZ ) *DOR 

C ROLL AND YAW OF THE MODEL 

IF(YZ.EQ.O. . AND.ZZ.EQ.O. ) THEN 
C CASE WHERE PITCH OF MODEL IS +/-90. 

PHIMT=ATAN2(-YX,ZX)*D0R 
PS1MT-0 . 

ELSE 

PHIMT-ATAN2(-YZ,ZZ)*D0R 
PS IMT-ATAN 2 ( -XY , XX ) *DOR 
END IF 
RETURN 
END 
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SUBROUTINE COMP (X,Y,Z,PSIOFF,THEOFF,PHIOFF,PSISB,THESB,PHISB, 
* THES ,PHIS ,XI ,YI ,ZI) 

C STING PITCH (Y) 

XB=X*COSD(THES)-Z*SIND(THES) 

YB*Y 

ZB=Z*COSD(THES)+X*SIND(THES) 

C STING ROLL (X) 

XC-XB 

YC”YB*COSD( PHIS )-ZB*SIND( PHIS ) 

ZC=ZB*COSD(PHIS)+YB*SIND(PHIS) 

C OFFSET YAW (Z) 

XD=XC*COSD(PSIOFF)-YC*SIND(PSIOFF) 

YD=YC*COSD( PSIOFF )+XC*SIND(PSIOFF ) 

ZD»ZC 

C OFFSET PITCH (Y) 

XE=XD*COSD(THEOFF)-ZD*SIND(THEOFF) 

YE=YD 

ZE=ZD*COSD(THEOFF)+XD*SIND(THEOFF) 

C OFFSET ROLL (X) 

XF=XE 

YF=*YE*COSD(PHIOFF)-ZE*SIND(PHIOFF) 

ZF=ZE*COSD(PHIOFF)+YE*SIND(PHIOFF) 

C STING BENDING IN YAW (Z) 

XG=XF*COSD(PSISB)-YF*SIND(PSISB) 

YG=YF*COSD(PSISB)+XF*SIND(PSISB) 

ZG=ZF 

C STING BENDING IN PITCH (Y) 

XH=XG*COSD(THESB)-ZG*SIND(THESB) 

YH=YG 

ZH=ZG*COSD(THESB)+XG*SIND(THESB) 

C STING BENDING IN ROLL (X) 

XI=XH 

YI=YH*COSD(PHISB)-ZH*SIND(PHISB) 

ZI=ZH*COSD(PHISB)+YH*SIND(PHISB) 

RETURN 

END 
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SUBROUTINE ALPBET ( U , V , W , PSIM , THEM , PHIM , ALP , BET ) 
DOR-57.2957795 

C YAW (Z) 

UA-U *COSD(PSIM)-V *SIND(PSIM) 

VA-V *COSD(PSIM)+U *SIND(PSIM) 

WA-W 

C PITCH (Y) 

UB-UA*COSD(THEM)-WA*SIND(THEM) 

VB-VA 

WB=WA*COSD(THEM)+UA*SIND(THEM) 

C ROLL (X) 

UC-UB 

VC=VB*COSD(PHIM)-WB*SIND(PHIM) 

WC-WB*COSD(PHIM)+VB*SIND(PHIM) 

C ALPHA AND BETA 

I F ( WC . NE . 0 . . AND . UC . NE . 0 . ) THEN 
ALP-ATAN2 ( WC , UC ) *DOR 
ELSE 
ALP-O. 

END IF 

IF(VC . LT.-l . )VC— 1. 

IF (VC .GT . l.)VC= 1. 

BET-ASIN(-VC)*DOR 

RETURN 

END 


SUBROUTINE CONV (X,Y,XSAVE,YSAVE) 

Xl-XSAVE 
Yl-YSAVE 
X 2=X 
Y2=Y 

IF(X2.EQ.X1) GO TO 3 
IF(Y2.EQ.YI) GO TO 4 
XK“( Y2-Y1 ) /(X2-X1 ) 

X-X2-Y2/XK 

XA*(XL+X2)/2. 

DX=ABS(X2-X1 ) 

IF(ABS(X-XA).GT.3.5*DX) X=XA+3 . 5*SIGN( DX, (X-XA) ) 
GO TO IOO 

3 X=X2-Y2 
GO TO IOO 

4 X“(X2+XI)/2. 

100 XSAVE=X2 

YSAVE=*Y2 

RETURN 

END 
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